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ABSTRACT 

In this paper, we investigate the extent to which observations of molecular clouds can correctly 
identify and measure star-forming clumps. We produced a synthetic column density map and a 
synthetic spectral-line data cube from the simulated collapse of a 5000 Mq molecular cloud. By 
correlating the clumps found in the simulation to those found in the synthetic observations, clump 
masses derived from spectral-line data cubes were found to be quite close to the true physical properties 
of the clumps. We also find that the 'observed' clump mass function derived from the column density 
map is shifted by a factor of ~ 3 higher t han the true clurnp mass function, due to projection of 



low-density material along the line of sight. Alves et al. (20071) first proposed that a shift of a clump 
mass function to higher masses by a factor ot '6 can be attributed to a star formation efficiency of 
30%. Our results indicate that this finding may instead be due to an overestimate of clump masses 
determined from column density observations. 

Subject headings: ISM: clouds - ISM: structure - stars: formation - stars: luminosity function, mass 
function 



1. INTRODUCTION 

It is believed that most, if not all, s tars form in a 
cluster environment (Lada & Lada' 2003 ) within highly- 
compressed cold clouds of dense molecular hydrogen gas 
and dust. When observing star-forming regions, we look 
for evidence of small, compact regions called "clumps" - 
localized dense regions in giant molecular clouds. These 
clumps are sites of star formation, and many will con- 
tinue to collapse to form stars. Understanding star for- 
mation regions requires a detailed understanding of the 
properties and evolution of clumps. 

Recent studies (i Rcid et al. 2010; C urtis fc Richer| 2010[ 
Shetty et al.pOlOffr'ineda et al ||200yl [Smith et°ar| 



have shown that limitations ot current observing tech 
niques make it difficult to correctly identify and mea- 
sure properties of these clumps and cores that reflect the 
true nature of the star-forming regions. Observations 
of star-forming regions in the night sky lack one funda- 
mental component: the third spatial dimension. It is 
important to understand how the projection of a three- 
dimensional structure can introduce bias into the anal- 
ysis of observations and affect subsequent conclusions. 
Issues which arise from projection effects, assumptions 
regarding cloud properties, and methods used to extract 
and analyse data can impact the outcome and interpreta- 
tion of results. While the Hcrschel and ALMA observa- 
tions will disentangle some of these issues, the most self- 
consistent way to understand the observational biases is 
by using large-scale simulations to model the collapse of 
a molecular cloud. 

In recent years, many groups have been in terested in 
studying star formation using simulations (e.g. Federrath 
et al.|2010l[Offner et al.|2009; Bate 2009; Tilleyfe Pudrltz 



2004p . These types ot simulations are able to provide a 



complete picture of the collapse of the molecular cloud, 
revealing the various effects of turbulence and gravity. 
Simulations of star-forming clouds can also provide in- 
sight into physical parameters without confusion from 



projection effects. 
Using the Sm oothed Particle H ydrodynamics (SPH) 



code. Gasoline, jPetitclerc] ( |2009[ ) modelled the collapse 
of a large 5000 M© molecular cloud. In order to make a 
direct comparison with observations, we produced syn- 
thetic observations from his results. We compare the 
clumps found using a common observer's tool - the 
publicly-available clum p-finding algorithm, Clumpfind 
(Williams et al. 1994) - on the three-dimensional sim- 
ulation to those found from the simulated observations 
to show the discrepancy between the properties of 'ob- 
served' clumps and 'actual' clumps. 

Clumpfind is one of many algorithms currently used 
for source extraction. Several new techniques have 
been developed in recent years, including getsources 
( Men'shchikov e t al. 2010, 2012) an d dendrogram meth- 
ods ( e.g. iRosoTowsky et al.||2068{ [Shetty et a"L]|2010 



Kauffmann et al. 2010p . An exploration of these alter 
nate methods ot clump identification is beyond the scope 
of this paper and has been left for future work. 

Studying the earliest stages of star formation can pro- 
vide clues as to the origin and universality of the stellar 
initial mass function (IMF), the distribution of stellar 
masses at the time of their birth. Since it has been shown 



that stars form from dense cores (e.g. 


Kirk et al. 


2006 


Motte et al. 


1998 


Johnstone et al. 


201 


) , many believe 



tribution of the clumps and cores. Motte et al. ( 1998 1 
extracted a clump mass function (CmI'') from observa- 
tions of the p Ophiuchi main cloud and were the first to 
no te the similar i ty to t he stellar IMF. In the seminal work 



by Alves et al 
IMl"' was a direc 



(|2007|), it was suggested that the stellar 
result of the CMF and that there exists 
a one-to-one correspondence between the two. This as- 
sumption has largely been adopted as the explanation for 
the resemblance of the two mass distributions. However, 



Reid et al. (2010) argue that the shape itself is a result 



of the central limit theorem and that this similarity is 
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Table 1 

Initial conditions of the SPH simulation 
of a giant molecular cloud 



Mass 


5000 Mq 


Particles 


36 088 472 


Diameter 


8 pc 


Initial uniform density 


300 cm-3 


Initial Joans mass 


2.8 Mq 


Particle mass 


1.3855 X 10""' 


Free-fall time 


1.9 Myr 


Temperature 


10 K 


Minimum Jeans mass 


0.01 Mq 


Minimum smoothing length 


50 AU 


Minimum spatial resolution 


100 AU 



M,7 



not physically significant. Alves et al. (20071 also argue 
that the two mass functions ditter only by a uniform star 
formation efficiency factor of approximately 30%. This 
is implied by a shift of the CMF to higher masses by a 
factor of 3. Although many are now questioning whether 
the IMF is universal across all star-forming regions and 
if the CMF h as a direct one-to-one correspondence with 
the IM F (e.g. "Smith et al.' 2008" 'Swi ft fc Wilf imls"2008'; 
Reid et al. 2010; Anathpindika 20 11, .Michel et al. 2011), 



no alternative explanation has been ottered tor the fac 



factor of 3 as a core-to-star conversion factor (e.f 




Enoch 


et al. 


2008 


Offner et al. 


2009 


Parmentier et al. 


2nn). 


in 


ihis paper, we wil 


explore these issues using our 



simulated observations. 



In section 2, we provide details 
simulation and in section 3, we 



of the ,Petitclerc, ( ,2009J 
explain the methods used to create our synthetic obser- 
vations from his results. In section 4, we present our 
results followed by a discussion of these results in section 
5. 

2. THE SIMULATION 
Our group has previously simula ted star formation in 



molecular clouds ( Petitclerc [2009 ) using the sm oothed 
particle h ydrodynamics (SPH) code. Gasoline (jWads 



ley et al.|[2004j . We modelled the collapse of an initially 
spherical 5000 Mq giant molecular cloud, initially 8 pc 
across, using approximately 36 million particles. The key 
specifications of the simulatio n are listed in Ta ble fT] We 
use an equation of state from Bate fc Bonnell| ( |2005| ) to 
model the thermal behaviour ol the gas. We have chosen 
this equation of state, as it mimics the behaviour of tem- 
perature and density seen in spherically-symmetri c col- 
lapses of m olecular clouds using radiative transfer (Bate 
et al. 2003). The value of lO'^^ 

g cm ^ in this equation 
ot state is taken to be the opacity limit f or fragmentation 
( [Bate et al.||2003[ [Bate fc Bonnell|2005| . Since our max- 
imum density of 10~^^ g cm~'^ is well below this density 
limit of 10~^^ g cm~'^, the entire simulation is isothermal 
at lOK and the gas is optically thin. In this paper, we 
present the results at 1.3 million years after the simu- 
lation began - approximately 70% of the initial free-fall 
time. 

These simulations had open boundary conditions, al- 
lowing us to simulate an entire molecular cloud undergo- 
ing gravitational collapse, rather than just a small por- 
tion within it. Other simulations have tended to focus 
solely on sub-regions within the clouds with masses only 



reaching as high as 500 solar masses (e.g. Bate 2009 
Smith et al.|2008[lTilley fc Pudritz|2004[ ). However, tEese 
sizes are typicarfor a si ngle clump seen in observations 



( [Bergin fc Tafalla 
rateiy represent w. 



2007 ). Thus, in order to more accu- 
lat IS seen in observations, we use sim- 



ulations of star formation that cover much larger scales. 
3. SIMULATED OBSERVATIONS 

In order to compare the results of the simulation with 
observations, we have created synthetic observations of 
a three-dimensional (position-position- velocity or 'PPV') 
spectral-line data cube and a two-dimensional (position- 
position or 'PP') column density map. These data are 
a more realistic representation of observations, as ob- 
servers typically do not know the physical position of an 
object along the line of sight. Our simulation provides 
the third spatial dimension necessary for comparison to 
the synthetic observations (position-position-position, or 
'PPP'). 

3.1. PPP 

To compare real objects found in our simulations to 
objects seen in our simulated observations, we need to 
characterize the simulation in a way that is complemen- 
tary to the observations. Therefore, we created a 3D 
smoothed density map of the gas in the simulation by 
interpolating the SPH particles to a 500 x 500 x 500 grid. 
The full spatial extent of the simulation is 3 x 10^ AU on 
a side so each pixel corresponds to 6000 AU. Any physi- 
cal quantity associated with the gas can be estimated by 
interpolating between the function values at the known 
positions of the pa rticles. A kernel function (Monaghan 
fc Lattanzio|1985 ) converts the particle information into 
a smoothed estimate. The weighting determines what 
fraction of the particle mass will contribute to each grid 
cell. The sum of the kernel weights is normalized to en- 
sure the total mass of the particle is counted only once. 
If the size of the particle carries it off the edge of the 
grid, it is assigned to the edge of the grid and if there 
are no grid lines intersecting the smoothing sphere, the 
particle weight will be assigned to the nearest grid cell. 
The cube can be rotated and "observed" along any line 
of sight. 

3.2. PP 

Star-forming regions are often observed in dust con- 
tinuum observations, which provides a measure of the 
column density of gas along the line of sight. Our two- 
dimensional (PP) column density map was produced 
from the simulation by integrating the density along the 
line of sight using the collapse task from the kappa soft- 
warepackage, a part of the Starlink software environ- 
menirl We arbitrarily chose the z axis of the simulation 
as our line of sight. The two-dimensional column den- 
sity map produced can be seen in Figure [T] We assume 
that the column density of material is directly analogous 
to the observed flux in sub-mm or IR observations since 
our gas is optically thin throughout the simulation time. 
Any heating which occurs in these relatively unevolved 
cores due to accretion is immediately radiated away. 

The spatial resolution is the same as the PPP cube: 
6000 AU/pixel and an image size of 500 x 500 pixels. 

^ http:/ /starlink.jach. hawaii.edu 
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Figure 1. Two-dimensional column density map of the simulated 
molecular cloud. The resolution of the map is 6000 AU/pixel. The 
total size of this box is 1.25 X lO" AU by 1.25 X lO** AU. Note 
the filamentary structure comparable to that seen in observations. 
The column density is in units of solar mass per pixel. 



The resolution was chosen to be comparable to that of an 
observation of the active Orion star-forming region using 
the SCUBA-2 instrument on the James Clerk Maxwell 
Telescope (JCMT). For the synthetic column density 
map with this resolution, the densest gas would corre- 
spond to a column density of approximately 1 Mq /pixel 
or lO'^"^ cm^^. These peak values are comparable to those 
seen in real observations, as starless cores have been ob- 
served to have typi cal column densities < 10^ 
( jCaselh et a"L]|2002a |. 



km s^^ per bin, which corresponds to a 3.050 kHz chan- 
nel separation for the NH3 spectrum. This was chosen 
to approximately ma tch that in the NH3 observations of 
a dense core used by Pineda et al. (2010). 

Since actual molecular lines have a tinite width, we 
must introduce a method to broaden our lines. We 
have chosen to use thermal broadening. To determine 
the width of our lines, we used the temperature of the 
isothermal gas in our simulation (10 K) and the molecu- 
lar mass of NH3. Using our velocity bin width and cor- 
responding frequency channel separation, as well as the 
rest frequency of NH3, we determined the width to be 7.6 
kHz or 0.1 km s~^. This width, together with the line- 
of-sight velocities (and corresponding frequencies) of the 
particles, allow us to broaden each of the spectral lines 
in our PPV cube. 

We produced a PPV cube where only the gas emis- 
sion traced by NH3 is included. The NH3 transition is 
excited in cl oud cores at a critical excitation density of 
1800 cm-3 ( |Tielens||2005| and NH3 is expect ed to de- 
plete only at very high densities > 10^ cm~^ (Bergin k, 
iLangcr 1_997'). NH3 is also abundant in dense cores and 
its continued presence at the high densities of core interi- 
ors makes NH3 is one of the best tracers of the phys ical, 
chemical, and k inematic properties of dense cores (Pa- 
gani et al.|2005 ). We used a threshold of 1800 cm~^ and 



did not include the contribution of any gas below this 
density when making this PPV cube. 

4. CLUMP-FINDING AND CORRELATIONS 

4.1. Studying Clumps in the Third Dimension 

With our simulation of molecular cloud collapse, we 
have 'observed' the clumps formed within t he cloud us- 
ing the c lump- finding algorithm, Clumpfind ( Williams et 



al. 



1994|), which we will refer to as Clumpfind3U when 
The Clumpfind algorithm requires 
which in observations is 



applied to 3D data, 
that we define a value for cr, 
the noise level. The algorithm then determines intensity 
contours in multiples of a for the entire dataset. The 
user specifies step increments and the minimum contour 
level in multiples of tr. The algorithm moves through the 
dataset in steps of the specified increment and assigns 
pixels to clumps based on the landscape of the local con- 
tour levels. The algorithm stops assigning pixels when it 
reaches the specified lowest contour level. 
The largest single value in a pixel for the PPP data 



cube was 0.8 Mj; 



To obtain an estimate of the 



3.3. PPV 

Submillimetre observations using spectral imag ing sys- 
tems, such as the HARP and ACSIS instruments ( Buckle 
et al. 2009) on the JCMT, allow astronomers to record 
a third dimension in the star-forming region. This third 
dimension is not a spatial dimension, but rather the ve- 
locity of the gas along the line of sight. This observa- 
tional technique produces images at various wavelengths 
to form a spectral-line data cube. 

To create the synthetic spectral-line data cube, we use 
the velocities of the particles to determine the motion of 
the gas along the line of sight (which we again take to be 
the z axis). For each particle, a contribution was added 
to the appropriate velocity bins (or channels) for all lines 
of sight the particle intersects. This produces infinitely 
sharp 'spectral lines'. For a given velocity channel, each 
contributing particle mass in that bin is divided up into 
grid cells by their specific weights. Combining each frame 
of column density per velocity channel results in a cube 
with two spatial dimensions and one velocity dimension. 

Our aim is to model molecular line emission along the 
line of sight; therefore, we require a molecular tracer 
to simulate the spectral-line emission. We have chosen 
the NH3 (1,1) hyperfine transition as our molecular line 
tracer. NH3 is an excellent trace r of dense cores, as it 

is present at very high densities (Caselli et al. 2002a|b 

Bergin &: Tafalla 2007). Our PPV cube is highly re- clumps are basically independent of noise levels, when 
solved with 148 velocity bins with a resolution of 0.04 the dynamic range of the observations are comparable to 



Clumpfind3D threshold for the PPP cube, we divided 
this maximum pixel value by a typical dynamic range. 
The dynamic range is defined for observations as the ra- 
tio between the brightest intensity and the background 
noise. As an example, the dynamic range on the Herschel 
instrument s, PACS and SPIRE , are ~1000 and -200, re- 
spectively ( Motte et al.| 2010|. A conservative estimate 
of 160 was used as the dynamic range for the PPP cube, 
resulting in a cr value of approximately 0.005 Mq/pIxcI 
for a resolution of 6000 AU/pixel. We used this dynam- 
ical range to set the Clumpfind threshold rather than 
adding noise to our simulations because the choices for 
the noi se properties wou ld be arbitrary. We have also 



shown (Reid et al. 2010) that the derived properties of 
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Table 2 

Clumpfind2d Parameter Sensitivity 



1000 



a 

(Mo/px) 


Low 


Step 
{^) 


# of clumps 


Mass in clumps 
(Mo) 


Percentage of 
Total Mass 


0.005 


So- 


2o 


886 


4937 


99% 


0.006 


So- 


2o 


558 


4777 


96% 


0.007 


So 


2o 


532 


4723 


94% 


0.006 


5o 


So 


455 


4558 


91% 


0.005 


5o 


So 


50S 


4648 


93% 



100 



that of modern instruments. 

Setting the lo-west density contour level in 
ClumpfindSD to 3(7 and the contour increment to 
2cr results in the identification of 499 clumps in the PPP 
data cube. These clumps contain approximately 430 
of material or 9% of the total mass in the simulation. 



4.2. Clump Correlations, Round One: PP vs. PPP 

We ran the t-wo-dimensional implementation of 
Clumpfind (-which we -will refer to as Clumpfind2D in 
this paper) on the synthetic column density map. All 
clumps -were identified using the findclumps task from 
the source extraction soft-ware package CUPID, -which is 
a part of the Starlink soft-ware collection. We used a va- 
riety of different parameters to determine their effect on 
the number of clumps identified by the algorithm. Our 
results are sho-wn in Table [2] -where 'Lo-w' is the lo-west 
contour level and 'Step' is the contour interval. Although 
changing the parameters can cause the number of iden- 
tified clumps to vary by a factor of up to 2, the amount 
of mass in clumps is independent of the parameters cho- 
sen. The mass in clumps makes up > 90 % of the mass 
of the entire cloud. Clumpfind2D decomposes the cloud 
into clumps by assigning almost all material to clumps. 
For this reason, many of these clumps may not necessar- 
ily connect to any physical structures -within the cloud. 
The inclusion of all emission along the line of sight, even 
lo-w-density gas unassociated -with the dense region, could 
affect the determination of clump properties. 

It has previously been sho-wn that clump pro perties 
determined from Cl umpfind2D can be unreliable (Kain- 
ulainen et al. 2009). Our results confirm that a small 
change of parameters in a clump-finding algorithm can 
alter both the number of clumps located and their prop- 
erties. It is difficult to obtain sets of parameters for 
the Clumpfind algorithm for the PP and PPP cases 
that identify the same clumps in both 'observations'. 
Rather than identify clumps independently in both PP 
and PPP, -we have taken the clumps identified in the 
three-dimensional volume and determined their position 
and area in the xy plane. We used a pixel mask made 
from this projection of clumps from the PPP cube and 
identified them directly onto the PP map. The mass 
of each PP clump includes all the mass along the lines 
of sight -within the area traced out by the PPP clump. 
This correlation method ensures that the clumps have 
an equal projected area and pixels -with identical (x,y) 
coordinates. 

Figure[2]compares the CMFs of clumps identified in the 
PP column density map using the Clumpfind2D method 
and our correlation method. The correlation method 
finds the projected masses of the PPP clumps by in- 
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Figure 2. The shaded region shows the range of possible mass 
distributions for clumps identified from the PP column density map 
using Clumpfind2D for the cases listed in Table 2. This is compared 
to a clump mass function of clumps identified from the same PP 
column density map using our correlation method (dashed line). 



eluding all the emission along the line of sight. The 
shaded region sho-ws the range of distributions produced 
by Clumpfind2D for the cases listed in Table [2] The 
range is tightly constrained at the high-mass end and 
the CMF produced by the correlation method provides 
a lo-wer limit to the range of possible mass distributions 
produced by the Clumpfind2D method. 

4.2.1. Global Properttes 

Clumpfind3D identified 499 clumps in the PPP cube 
and returned the integrated intensity and the peak (x,y,z) 
positions of each of these clumps. The total mass con- 
tained in clumps in the PPP cube is approximately 430 
M0. Using the correlation method, the total mass con- 
tained in clumps from the masked PP map, -which has 
an equal area in the XY-plane to the clumps in the PPP 
cube, is 1647 M0. This is not surprising as large amounts 
of low density material along the line of sight are no-w 
contributing to the mass estimate of the clumps in those 
columns (see Figure [s]) . When the simulation is inte- 
grated along this axis to produce a t-wo-dimensional col- 
umn density map, all of this emission contributes, result- 
ing in an identified object -with a mass -which is larger 
than the real star-forming clump. In some cases (e.g. 
the bottom panel of Figure pi) there is no clump identi- 
fied in the PPP cube, but tne integrated intensity of a 
significant amount of low-density gas can masquerade as 
a clump in the projection on the plane of the sky. 

Figure |4] shows the mass distribution of the 499 clumps 
for both the PPP case (solid line) and the PP projection 
(dashed line). Both clump mass functions closely reflect 
the power-law t rend at high mass seen in the stellar IMF 
( Sal£eter| 1955), shown by the dash-dotted line i n Fig- 
We have also plotted the Alves et al. ( 2007 ) dense 



core mass function for the Pipe JNebula, shown m Fig- 
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Figure 3. Cross-section slices of the simulation in the xz-plane. The colour scale shows density in units of cm~^. Large filaments of 
low density gas are present along the z-axis, chosen to be our line of sight, when observed in the xz-plane. The top and middle panels 
respectively show the largest clump and a typical clump (both circled) identified in the PPP cube plus the foreground and background 
material which contributes to their PP masses. These clump masses are shown as circled points in Figure [5] The bottom panel shows a 
region which was identified as a clump in two dimensions using Clumpfind2D, but contains no true clump as identified by ClumpfindSD. 
The white dashed lines contain the lines of sight which intersect the circled clump of interest, where the spacing of the white lines is 
determined by the size of the clump. These slices have a width of 1. 5 X 10^ AU and a length of 1.5 X 10^ AIJ, extending fully across the 
simulation. This visualization was made using SPLASH |jPriH^(2007)l. 

urelHas circles with error bars. We note that the factor 
of « 3 difference in mass between the PP and PPP clump 
masses is due entirely to the contribution of low-density 
gas and multiple cores along the line of sight in the PP 
case. While our PPP clump mass function is well-fit by 
a Salpeter mass function, these objects are still clumps, 
not protostars. The mass difference here cannot be in- 
terpreted as a star formation efficiency. 



1000 



4.2.2. Clump-by- Clump Comparison 

Using the masses obtained for each of the 499 clumps 
from our correlation method, we have plotted the ob- 
served masses of each clump (in PP) against their cor- 
responding actual masses (in PPP) to produce a direct 
clump-to-clump comparison. The power-law trend ob- 
served between Mppp and Mpp in Figure [5] exists due 
to the correlation between Mppp and the projected area 
of the PPP clumps. Though the column of material along 
the line of sight (as shown in Figure^ is not correlated 
with the PPP masses, as can be seen from the scatter in 
Figure [5] when multiplied by the area, a net correlation 
arises. 

The functional form 

Mpp = (6 ± 1) (Afppp)°-^*°-' 

was fit to the data with a reduced chi-squared of 3.5 x 
10~^. The slope of this relationship is w 2/3, and we ex- 
pect this is related to the ratio of surface area (cx R^) to 
volume ((X R'^) of the clump, since the density is approx- 
imately constant. For a given PPP mass between 2M0 
and IOM0, this function can provide a range of possi- 
ble PP masses that are approximately 3-5 times greater. 
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Figure 4. Clump mass function of the clumps identified in the 
PPP data cube (solid line). Also shown is the clump mass function 
of the same clumps with their masses estimated from the 2D pro- 
jected PP image using the correlation method (dashed line). The 
dash-dotted line is the Salpeter (1955) IMF. The dense core mass 
function (circles) obtained by ^Alves et al.^ ( ,2007^ from observations 
of the Pipe Nebula is also shown tor comparison. 



which is again comparable to the factor of three often 
quoted in other works with regard to the CMF ( Alves 
et al.||2007l [Enoch et alTjpOOSl |Offner et aL]|2009j 'I'he 
estimated PP masses can be sven more than 3-5 times 
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Figure 5. Observed mass from the column density map vs. the 
actual mass from the simulation cube. A fit to the data is shown 
as a solid line, and given on the figure. The dashed line shows the 
expected correlation if the observed mass were not contaminated 
by foreground and background material. The circled points are the 
clumps shown in the top and middle panels of Figure |3] 

larger than the actual mass for lower mass objects. 

4.3. Clump Correlations, Round Two: PPV vs. PPP 

A second method of determining properties of star- 
forming clumps is to observe them using spectral line 
data. We follow the same method of identifying clumps 
in the PPP cube and then studying the properties of the 
same region in the PPV 'observation'. In order to ob- 
tain a PPV mask of the desired clumps from the PPP 
cube, we first obtained the (x,y) pixel coordinates of the 
PPP clumps to create a 2D mask. We then determined 
the mean velocities of the PPP clumps to correlate the 
clumps in (x,y,z) coordinate space with (x,y,v^) coordi- 
nate space. 

The mean velocity and velocity dispersion of each of 
the PPP clumps were determined using the velocities of 
the particles in the simulation. We are only concerned 
with the mean velocity along the line-of-sight, v^, for 
comparison with the PPV cube. Once the mean veloc- 
ity and velocity dispersion of each clump are known, 
we use these values to estimate a range of velocities, 
[vz — az,Vz + (Tz], to represent the range of velocity chan- 
nels spanned by a Doppler-broadened line profile in the 
velocity spectrum. 

For each clump, we specify the velocity bins we would 
like to extract from the PPV cube. We also specify 
the minimum volume density of gas that can contribute 
to the NH3 emission (1800 cm~'^). The mass of the 
PPV clumps is the contribution of all particles in high- 
density regions moving at line-of-sight velocities between 
Vz — and Vz + Uzi with the same xy area as the PPP 
clumps. We have tested this method of clump selection 
against the clumps found by applying ClumpfindSD di- 
rectly on the PPV cube, and the clump mass functions 
are very similar. Therefore, we analyze the properties of 
the clumps found by our correlation method in the PPV 
cube. 

4.3.1. Global Properties 

The total mass contained in clumps in the PPP cube 
is approximately 430 M©. Using the correlation method 
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Figure 6. Clump mass function of the clumps identified in the 
PPP data cube (solid line) compared to the clump mass function of 
the same clumps with their masses estimated from the NH3 (1,1) 
emission of the {x,y,V2) PPV data cube (dashed line). 



to locate these clumps in the PPV data cube, we find 
that the total mass contained in clumps from the masked 
PPV cube is 439 M©. We have produced a clump mass 
function for the 499 clumps identified in the NH3 PPV 
data cube, shown in Figure |6] 

The PPV clump mass function is extremely close to 
the true PPP clump mass function. We performed a 
two-sided Kolmogorov-Smirnov Test (KS-Test) to com- 
pare the PPV and PPP masses. There was an 81% prob- 
ability that the PPV and PPP masses were drawn from 
the same distribution. Most of the extended emission, 
even that moving at a similar line-of-sight velocity as 
the clump, does not contribute significantly in this case. 
Based on these results, it appears that clump properties 
derived from PPV spectral-line data cubes are quite rep- 
resentative of the true physical properties of the clumps. 

4.3.2. Clump-by- Clump Comparison 

The comparison of actual and observed masses using 
the PPV spectral cube are shown in Figure [7] Using 
the masses obtained for each of the 499 clumps from 
the correlation method, the observed PPV masses are, 
on average, the same as the actual PPP masses. We 
conclude that using spectral line data is the best method 
for accurate estimation of clump properties. 

5. DISCUSSION 

We have 'observed' a SPH simulation of the collapse of 
a spherical molecular cloud, using the techniques and 
tools common to real observations of star-forming re- 
gions. We determine the masses of three-dimensional 
star-forming clumps as they would be observed in a two- 
dimensional column density map and in a PPV spectral- 
line data cube. We found that projection has the great- 
est effect on estimates derived from PP column den- 
sity maps. We have seen that there are some objects 
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Figure 7. Observed mass (PPV) vs. the actual mass (PPP) 
of clumps. The PPV masses are obtained from only the highest 
density gas traced by NH3. The dashed line is the one-to-one line. 



identified in column density maps as clumps, which are 
not clumps at all. Instead, they are extended regions 
of low density gas which, when projected onto a two- 
dimensional sky, appear as a single object. We found 
that effects of projection can result in overestimates of 
clump masses derived from two-dimensional column den- 
sity maps. Our results show that the observed shift of 



the CMF f rom the stellar IMF to higher masses (Alyes 
et al.|2007 ) is not the result of a star formation efficiency 
factor. 

By correlating the clumps found in PPP to those found 
in PPV, we find that the properties of objects derived 
from PPV spectral-line data cubes were more represen- 
tative of the true physical properties of the clumps than 
those obtained from the PP column density map. We 
conclude that clump properties are best derived from 
molecular-line data rather than continuum data, in order 
to minimize projection effects and maximize the identi- 
fication of truly bound structures in observations. 

In addition to the work presented here, we investigated 
various modifications of our simulation with different 
evolutionary stages, different orientations, and for dif- 
ferent initial collapse conditions. Recent studies (|Arzou- 
manian et al.||20lT| [Schmalzl et~ar]|20T0l |lVLen'shchikov 
et al.||2010p have shown that filaments are very common 
m molecular clouds and may have a significant role in 
star formation. To learn more about the importance of 
initial geometry, additional simulations were run with 
different large scale turbulent velocity modes included in 
the initial conditions. When we study the effect of the 
initial turbulent velocity spectrum on the formation of 
dense cores, the anisotropic collapse cases all found sig- 
nificantly more dense clumps and more mass contained 
in clumps than that seen for the spherical collapse case. 
This result implies that dense cores in the cloud are more 
inclined to form in filamentary-like structure, which is in 
agreement with conclusions d rawn from recent observa- 
tions of star- forming regions ([Arzoumanian et al. 2011 



Schmalzl et~aL]|2010[ |Men'shchikov et al.||2UlUp . How- 



are unchanged, regardless of the evolutionary stage, the 
orientation, and the initial collapse condition. 

The next step in this work is to look at an even larger 
(and more realistic) simulation. We are currently run- 
ning a simulation of a 50000 Mq molecular cloud to 
track the evolution of the cores from the initial collapse of 
the cloud through to the pre-main sequence evolutionary 
phase. As the cores evolve, we will be able to see how of- 
ten they will merge and collapse into bound objects and 
how often they disperse to determine a statistical like- 
lihood that an observed core will eventually become a 
star. The study of the core evolution will provide insight 
into the formation mechanisms of multiple star systems, 
such as binaries, and will also give an indication of the 
multiplicity of these objects in observations; for exam- 
ple, how often a single, identified core is, in reality, two 
or more widely separated cores which have been confused 
in projection. 

We also intend to run a simulation incorporating the 
effect of radiative feedback from the protostars using a 
radiat ive transfer code recentl y implemented in Gaso- 
line (Rogers & Wadsley 2011). Since we are interested 
in studying the later stages of pre-main sequence stellar 
evolution when the cores become associated with a com- 
pact source of luminosity and are emitting protostellar 
radiation, we need to understand the effects of localized 
heating due to the radiative feedback from newly-formed 
stars and determine the significance of the addition of ra- 
diative transfer under the conditions of our simulation. 

ACKNOWLEDGMENTS 

We would like to thank SHARCNET (Shared Hier- 
archical Academic Research Computing Network) and 
Compute/Calcul Canada, which provided dedicated 
resources to run these simulations. We would also like 
to thank Christine Wilson and Michael Reid for their 
helpful conversations. This work was supported by 
NSERC. JW acknowledges support from the Ontario 
Early Researcher Award (ERA). 



REFERENCES 



Alves, J., Lombard!, M. & Lada, C. J. 2007, A&A, 462, L17 
Anathpindika, S. 2011, New Astronomy, 16, 477 
Arzoumanian, D., Andre, P., Didelon, P., et al. 2011, A&A, 529, 
6 

Bate, M. R. 2009, MNRAS, 392, 590 
Bate, M. R. & Bonnell, I. A. 2005, MNRAS, 356, 1201 
Bate, M. R., Bonnell, I. A., & Bromm, V. 2003, MNRAS, 339, 577 
Bergin, E. A. & Tafalla, M. 2007, ARA&A, 45, 339 
A. & Langer, W. D. 1997, ApJ, 486, 316 

Hills, R. E., Smith, H., et al. 2009, MNRAS, 399, 



Bergin, E 
Buckle, J. V 

1026 
Caselli, P. 

331 
Caselli, P. 



Walmsley, C. M., Zucconi, A., et al. 2002a, ApJ, 565, 
Myers, P. C, 



ever, we found that our results regarding the relation- 
ship between the 'observed' masses and 'actual' masses 



Benson, P. J., Myers, P. C, & Tafalla, M. 2002b 
ApJ, 572, 238 
Curtis, E. I. & Richer, J. S. 2010, MNRAS, 402, 603 
Enoch, M. L., Evans, N. J., Sargent, A. I. et al. 2008, ApJ, 684, 
1240 

Federrath, C, Banerjee, R., Clark, P. C, & Klessen, R. S. 2010, 
ApJ, 713, 269 



8 



Ward et al. 



Johnstone, D., Wilson, C. D., Moriarty-Schicven, G., et al. 2000, 
ApJ, 545, 327 

Kainulainen, J., Lada, C. J., Ratliborne, J. M., & Alves, J. F. 

2009, A&A, 497, 399 
Kauffmann, J., Pillai, T., Shetty, R., Myers, P. C, & Goodman, 

A. A. 2010, ApJ, 712, 1137 
Kirk, H., Johnstone, D., & Di Francesco, J. 2006, ApJ, 646, 1009 
Lada, C. J. & Lada, E. A. 2003, ARA&A, 41, 57 
Men'shchikov, A., Andre, P., Didelon, P., et al. 2010, A&A, 518, 

103 

Men'shchikov, A., Andre, P., Didelon, P., et al. 2012, 

arXiv: 1204.4508 
Michel, M., Kirk, H., & Myers, P. C. 2011, ApJ, 735, 51 
Monaghan, J. J. and Lattanzio, J. C. 1985, A&A, 149, 135 
Motte, F., Andre, P., & Neri, R. 1998, A&A, 336, 150 
Motte, F., Zavagno, A., Bontemps, S., et al. 2010, A&A, 518, L77 
Oflfner, S. S. R., Klein, R. I., McKee, C. F., & Krumholz, M. R. 

2009, ApJ, 703, 131 
Pagani, L., Pardo, J.-R., Apponi, A. J., Bacmann, A., & Cabrit, 

S. 2005, A&A, 429, 181 
Parmentier, G., Kauffmann, J., Pillai, T., & Menten, K. M. 2011, 

MNRAS, 416, 783 
Petitclerc, N. 2009, PhD Thesis, McMaster University 



Pineda, J. E., Goodman, A. A., Arce, H. G., et al. 2010, ApJ, 
712, L116 

Pineda, J. E., Rosolowsky, E. W., & Goodman, A. A. 2009, 

ApJL, 699, L134 
Price, D. J. 2007, PASA, 24, 159 

Reid, M. A., Wadsley, J., Petitclerc, N., & Sills, A. 2010, ApJ, 
719, 561 

Rogers, P. D. & Wadsley, J. 2011, MNRAS, 414, 913 
Rosolowsky, E. W., Pineda, J. E., Kauffmann, J., & Goodman, 

A. A. 2008, ApJ, 679, 1338 
Salpeter, E. E. 1955, ApJ, 121, 161 

Schmalzl, M., Kainulainen, J., Quanz, S. P., et al. 2010, ApJ, 
725, 1327 

Shetty, R., CoUins, D. C, Kauffmann, J., et al. 2010, ApJ, 712, 
1049 

Smith, R. J., Clark, P. C, & Bonnell, I. A. 2008, MNRAS, 391, 
1091 

Swift, J. J. & Williams, J. P. 2008, ApJ, 679, 552 

Tielens, A. G. G. M. 2005. The Physics and Chemistry of the 

Interstellar Medium, Cambridge University Press 
Tilley, D. A., & Pudritz, R. E. 2004, MNRAS, 353, 769 
Wadsley, J. W., Stadel, J., & Quinn, T. 2004, New Astronomy, 9, 

137 

WiUiams, J. P., do Geus, E. J., & Blitz, L. 1994, ApJ, 428, 693 



